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Interacting bosons with two “spin” states in a lattice show novel superfluid-insulator phase transitions in 
the presence of spin-orbit coupling. Depending on the parameter regime, bosons in the superfluid phase can 
condense to either a zero momentum state or to one or multiple states with finite momentum, leading to an 
unconventional superfluid phase. We study the response of such a system to modulation of the optical lattice 
potential. We show that the change in momentum distribution after lattice modulation shows distinct patterns in 
the Mott and the superfluid phase and these patterns can be used to detect these phases and the quantum phase 
transition between them. Further, the momentum resolved optical modulation spectroscopy can identify both the 
gapless (Goldstone) gapped amplitude (Higgs) mode of the superfluid phase and clearly distinguish between the 
superfluid phases with a zero momentum condensate and a twisted superfluid phase by looking at the location 
of these modes in the Brillouin zone. We discuss experiments which can test our theory. 


I. INTRODUCTION 

Ultracold atoms have emerged in recent years as a valuable 
platform to study strongly interacting many-particle Hamil¬ 
tonians relevant to condensed matter systems, nuclear matter 
and many other different fields of physic^iEl. The unprece¬ 
dented control over the Hamiltonian parameters and easy ac¬ 
cess to strongly interacting regimes have opened up the pos¬ 
sibilities of systematically studying interacting many body 
models which have been used as paradigms to describe a 
multitude of phenomena from the realms of condensed mat¬ 
ter physics, nuclear physics, astrophysics and high energy 
physics. In the arena of condensed matter physics, many body 
lattice Hamiltonians like the Bose- and the Fermi-Hubbard 
models, which are used as paradigms to study strongly in¬ 
teracting bosons and fermions, have been experimentally im¬ 
plemented and studied®. These provided a wealth of in¬ 
formation which is relevant to the phenomena of superfluid- 
insulator transitiorP and high temperature superconductor^ 
respectively. In addition, various spin-models have also been 
realized using ultracold atom^®, which will be useful in study 
of frustrated magnetic systems and spin liquid^. 

Gauge fields and their interaction with matter are at the 
heart of our understanding of the physical phenomena around 
us. While the electromagnetic interactions governing most 
of condensed matter physics is described by a U(l) Abelian 
gauge theory, more complicated non-Abelian versions gov¬ 
ern the weak and strong interactions. In condensed matter 
systems, the presence of externally imposed gauge configura¬ 
tions, e.g. a fixed electric or magnetic field, can lead to inter¬ 
esting and qualitative change in properties of the system, e.g. 
in presence of magnetic fields, vortex lattices can form and 
melt in a superconductoiE^, the ground state of the system can 
have non-trivial topology and associated quantized conduc¬ 
tance in quantum Hall effect etcp2. Non-Abeli an g auge fields, 
which can take the form of spin-orbit coupling^^is an essen¬ 
tial ingredient in the realization of topological insulator^ and 
topological superconductor^^! and plays a crucial role in un¬ 
derstanding novel phenomena like anomalous quantum Hall 
effect (AQHE) in systems with strong spin orbit coupling. 

Ultracold atoms have been dressed by laser fieldJi^, so that 


in the lowest manifold of dressed states, the effective Hamilto¬ 
nian is that of the bosons/fermions, whose spin a nd orb ital de¬ 
grees of freedom are coupled. Alternate proposal^i^E^of real¬ 
izing effective spin-orbit coupling terms are present in the lit¬ 
erature. Recently time-varying magnetic field gradients have 
been used to realize spin-orbit couplings! The implementa¬ 
tion of specific types of spin-orbit coupling can lead to inter¬ 
esting phases of matter like topological insulators and topo¬ 
logical superconductors. The ability to tune both the effective 
spin-orbit coupling and the interaction strength in these sys¬ 
tems, which is very hard to achieve in material based systems, 
has opened up the possibility of studying novel Mott insu¬ 
lators and superfluids with Bose Einstein condensation into 
states with finite momenta in these system^S®. 

Compared to the wide array of experimental techniques 
available to probe material systems, cold atom systems suf¬ 
fer from a paucity of experimental probes. The main tool 
of obtaining time of flight absorption images, which trans¬ 
lates to observation of momentum distribution in lattice sys¬ 
tems, is a rather blunt instrument to differentiate between the 
myriad phases of matter which can occur in these systems. 
Eurther, the simple time of flight measurement does not pro¬ 
vide any dynamic (energy dependent) information about the 
system, which is crucial in understanding its low tempera¬ 
ture properties. A few spectroscopic techniques like rf spec- 
troscop v^, B ragg spectroscopjE^and lattice modulation spec- 
troscopjEUD are available to obtain energy resolved informa¬ 
tion about these systems. The latter method constitutes an 
ultracold atom counterpart of standard angle-resolved photoe¬ 
mission spectroscopy and provides energy and momentum re¬ 
solved information regarding the single particle spectral func¬ 
tion of the bosons. This method has been proposed for single 
species bosons in the strong-coupling regim^^, however to 
the best of our knowledge, it has never been applied to spinor 
bosonic systems with spin-orbit coupling. Here, we will focus 
on the lattice modulation spectroscopy of such systems. 

In this paper, we will consider two component boson^^ZRS] 
in a 2D square optical lattice, which are interacting with a lo¬ 
cal Hubbard type interaction. We will consider “spin”-orbit 
couplinin these systems, implemented either by Raman 
dressing of the atoms or modulation of magnetic held gradi- 
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ents. We consider the response of this system to optical lat¬ 
tice modulation spectroscopy and demonstrate the following. 
First, we show that optical lattice modulation spectroscopy 
can resolve the Mott and the superfluid phases both by look¬ 
ing at the absence/presence of gapless Goldstone modes and 
by identifying the unique pattern of appearance and disappear¬ 
ance of excitation contours in the Brillouin zone. Second, we 
demonstrate that the momentum resolved nature of the optical 
modulation spectroscopy can be used to clearly distinguish su¬ 
perfluids with condensate at zero momentum from those with 
condensate at finite momentum. Third, we provide a general 
theory of extracting the spectral function of spin-orbit coupled 
bosons in the strong coupling regime and near their superfluid- 
insulator phase transition by computing their response to the 
lattice modulation. Finally, we use the results of this theory 
to make definite predictions about the excitation spectrum of 
these bosons both in the Mott and the superfluid phases which 
can tested in realistic experiments and demonstrate that lattice 
modulation spectroscopy can reliably identify and character¬ 
ize both the gapless (Goldstone) and the gapped amplitude 
(Higgs) modes of spin-orbit coupled superfluid bosons; we 
note that such an analysis of the properties of excitation of 
the SF and Mott phases of these systems is beyond the scope 
of standard time-of-flight measurement which can also distin¬ 
guish between the Mott and the SF phases. 

The plan of the rest of the paper is as follows. In Sec.|I^ we 
discuss the general theory of optical modulation spectroscopy 
for spin-orbit coupled bosons. This is followed by Secs. |m] 
and 1^ where we inspect the detailed response of the system 
to optical modulations in the Mott and the superfluid phases 
respectively. Finally, we discuss possible experiments that can 
be carried out to validate out theory, sum up our main results, 
and conclude in Sec.lVl 




:: 
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FIG. 1. Density plots showing (a) intraband and (b) interband 
particle-hole excitation spectrum across the Brillouin zone in the 
Mott phase. The parameters of the system are J = 0.0217, 'y = 
0.0117, p = 0.217,17 = 0.0117 and = 0.4. The intraband ex¬ 
citations have energies between 0.2417 to 0.5717, while the inter¬ 
band excitations have energies between 0.7317 and 1.217. (c) Intra¬ 
band and (d) Interband particle-hole excitations in the Mott phase for 
J = 0.0117, 7 = 0.0417, p = 0.217,17 = 0.0117 and C = 0.4. In this 
case minimum of intraband excitations are shifted from the zone cen¬ 
ter to (Ico, fco) with feo = 1.131a“^. The figure shows the contours 
in the interval 0 < fca,, fcj, < tt; the contours in the other quadrants 
can be obtained by the reflection of the present figure around and 
ky axes and the origin. The intraband excitations have energies be¬ 
tween 0.2817 to 0.517, while the interband excitations have energies 
between 0.8717 and 1.117. Note that the color scheme has same ab¬ 
solute value in all the panels to show which excitations overlap with 
each other in energy. 


II. LATTICE MODULATION SPECTROSCOPY WITH 
NON-ABELIAN GAUGE EIELDS 


In a system with multiple species of bosons, Raman lasers 
can be used to generate non-abelian gauge fields which cou¬ 
ple the different species of bosons. In most experimental 
situations, the different boson species are actually different 
hyperfine state of the same bosonic species allowing one to 
treat them as a system of multi-component bosons. In fact, 
if one can trap two bosonic states it can be considered as an 
effective pseudo-spin 1/2 system, which is however made of 
bosonic atoms. A particularly interesting non-abelian gauge 
field configuration is the Rashba spin-orbit coupling in the 
two component bosonic system. There are various propos¬ 
als to implement the Rashba spin-orbit coupling, although 
currently experiments have focussed on the more easily im- 
plementable configuration of equal Rashba and Dresselhouse 
coupling, which leads to spins coupling to momenta in a par¬ 
ticular direction. Interacting bosons with Rashba spin-orbit 
coupling shows interesting chiral Mott and superf luid p hases 
as various parameters are tuned in the HamiltoniarPSSI 
The Hamiltonian for two-species bosons in a square optical 
lattice in the presence of Rashba spin-orbit coupling term can 


be written a} j 2Z| 28 | 34 | 

Uq — ^ ^ [( 17(7 1 )] QU 

i(T 

-J -f i7 ^ T-jz • (ct X 4) 4 (1) 

<U>o- (ij) 

Here bia annihilates a boson of spin a =t, i on the i*^ site, 
n-ia = bl^ha is the number of cr bosons, U{C,U) is the intra- 
(inter-)species interaction strength between the bosons, and 
J denotes the nearest neighbor hopping amplitude. Here p, 
and fl are the species independent and the species dependent 
chemical potentials; the latter acts as an effective Zeeman 
magnetic field for these bosons. The last term represents the 
lattice analogue of t he R ashba spin-orbit coupling generated 
by the Raman laser^SEl, vvith a coupling constant 7. Here, 
dij is unit vector along the x — y plane between the neigh¬ 
boring sites i and j, a is the vector of Pauli matrices, and 
= [bp, bii) is the two-component boson field. 
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(a) «= 0.36 U (b) co= 0.46 U 



FIG. 2. Lattice Modulation response in the Mott phase for dif¬ 
ferent modulation frequencies. The parameters of the system are 
J = 0.0217, 7 = 0.0117, = 0.217,G = 0.0117 and C = 0.4. 

As the frequency of modulation crosses the Mott gap, a contour of 
intrahand excitations are seen to disperse across the Brillouin zone ( 
a and h). Then, as the frequency crosses the interhand threshold, an¬ 
other contour of interband excitations disperse across the Brillouin 
zone (c and d). Note that the spectral weight of the interband transi¬ 
tions are much smaller than the weights of other transitions. 


The non-interacting part of the Hamiltonian is given by 

k 

A(fc) = efel + 7fcCT++7fc(T“ (2) 

where k = k = {kx,ky) is the 2 D quasi-momentum, = 
—2 J(cos kx+cos ky) and 7^ = —27i(sin kx — i sin ky). This 
can be diagonalized to obtain chiral bands touching each other 
at the zone center for H = 0, while a finite H opens up a gap 
between the two bands. The ratio 7/J controls the bare dis¬ 
persion including the location of the band-minima. As 7/J 
becomes larger the location of the minimum of the lowest 
band shifts away from the zone center, [0,0]. 

The spin-orbit coupled bosons undergo a Mott insulator to 
sup erfluid quantum phase transition as a function of J/U and 
7 /L ^^ 0 | 29 | Jjj jjjg strongly interacting limit, each site has ex¬ 
actly bosons of spin a, with the value of decided by /r 
and fl. The system is thus a Mott insulator with no number or 
spin fluctuations. As J/U or 7/17 increases, the system un¬ 
dergoes a phase transition to a state with delocalized bosons. 
This is the superfluid state with phase coherent Bose conden¬ 
sate and gapless Goldstone excitations due to broken 17 ( 1 ) 
symmetry. If the transition takes place when 7/J is large, 
the minimum of the effective dispersion occurs at [i/co, 
and the bosons condense into these finite momentum states, 
leading to a twisted superfluid phasd^. At low values of 7/J 
one recovers the standard superfluid with a condensate at zero 
momentum. Thus this system shows two remarkable qualita¬ 
tive changes : (a) a superfluid-Mott insulator transition as a 
function of J/U and 7/U and (b) a change from a standard 


(a) 0 )= 0.3 U 



FIG. 3. Lattice Modulation response in the Mott phase for dif¬ 
ferent modulation frequencies. The parameters of the system are 
J = O.OlU, 7 = 0.04U, p = 0.2U,O = O.OlU and C = 0.4. 
In this case minimum of intraband excitations are shifted from the 
zone center to (±lco, ±fco) with Icq = 1.131a“^. The figure shows 
the contours in the interval 0 < kx,ky < tt; the contours in the 
other quadrants can be obtained by the reflection of the present figure 
around k^ and ky axes and the origin. As the frequency of modula¬ 
tion crosses the Mott gap, a contour of intraband excitations are seen 
to disperse across the Brillouin zone starting from a circle around 
[kojko] ( a and b). Then, as the frequency crosses the interband 
threshold, another contour of interband excitations disperse across 
the Brillouin zone starting from the zone center (c and d). 


superfluid to a twisted superfluid as a function of 7/J. As 
we will show, lattice modulation spectroscopy can distinguish 
both these phenomena and hence provide us with a wealth of 
information about this system. 

The optical lattice modulation spectroscopy protocol that 
we propose consists of the following steps: (a) The optical 
lattice potential is weakly modulated with an a.c. field on top 
of the static field that forms the lattice in the original system, 
with the Raman fields and the trapping potential turned on. 
This leads to a modulation of the hopping parameter and the 
spin-orbit coupling, (b) The modulation is turned off after 
some time, making sure that the system is still in perturbative 
regime. At the same time, the Raman lasers, the optical lat¬ 
tice lasers and the trapping potential is also turned off and the 
system undergoes ballistic expansion, from which the (spin- 
resolved) momentum distribution of the system right after the 
modulation can be measured, (c) The change in the momen¬ 
tum distribution (from the unperturbed/ unmodulated system) 
will provide us with information about the spectrum and spec¬ 
tral weight of one particle excitations in these systemJ^. 

In experiments, the lattice Hamiltonian Eq. is imple¬ 
mented by putting a system of bosons, characterized by mass 
m, a continuum spin-orbit coupling 7c and an effective mag¬ 
netic field h generated by the Raman lasers, under a pe¬ 
riodic optical potential Vo[cos^(a:/a) + cos^(y/a)], where 
a = Aop/( 27 r) is a transverse confinement scale used to con- 
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(b) 6=71/4 



FIG. 4. Lattice modulation response in the Mott phase with the variation in frequency for various cuts in the Brillouin zone going radially 
outwards at an angle 6 with the axis. In (a) - (c) the parameters of the system are J = 0.0217, 7 = 0.0117, ^ = 0.217,fl = 0.0117 and 
^ = 0.4 with the band minimum at the zone center, (a) 6 = 7 r/ 2 , (b) 6 = 7 r /4 and (c) 6 = n/6. In (d) - (f) the parameters of the system are 
J = 0.011/, 7 = 0.041/, jj, — 0.21/,fl = 0.011/ and ( = 0.4. In this case minimum of intraband excitations are shifted from the zone center 
to (fco, fco) with ko = 1.131a“^. (d) 6 = n/2, (e) 8 — 7 r /4 and (f) 9 = 7 r/ 6 . 


Struct the 2D square lattice, and Aop is the wavelength of the 
laser forming the optical lattice. The interaction between these 
bosons in continuum are characterized by a intra-species scat¬ 
tering length Qs and an inter-species scattering length In 
the limit of a deep lattice, the continuum parameters are re¬ 
lated to the lattice parameters by 


U = 4^VoERa,/a, 

h 


(3) 


where Er = /2m is the recoil energy. It is evi¬ 

dent that J and 7 both depend exponentially on the lattice 
depth Vo, while U has only polynomial dependence on lattice 
height. Thus, an ac modulation put on the optical lattice depth 
Voit) = Vo + 5V cos (jjt will lead to simultaneous modulation 
of the hopping parameter, the spin orbit coupling as well as 
the interaction parameter U. Using U as an overall scale for 
the problem, the perturbation Hamiltonian, to linear order in 


variations of the optical lattice potential, can be written as 
H,=5UHo-U5(^\ 

(4) 

(G> 

where the first term, proportional to the variation of U, com¬ 
mutes with the unperturbed Hamiltonian and hence does not 
create excitations. This term can thus be neglected as far as 
modulation spectroscopy is concerned. In this case, the per¬ 
turbation consists of modulating the hopping and spin-orbit 
terms with amplitude^ 

X = SJ/J = 6j/j=^iSV/^/Vo^) (5) 

The perturbation Hamiltonian can then be written as 

Hi{t) = Acoswf E T'^^A(fc)T'fc. (6) 

k 


In linear response regime, the momentum distribution right 
after the modulation is turned off oscillates with the frequency 
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of the perturbatioiJSl 

Sna,k (t) = 6n^^\a,k,uj) cos ujt + {a,k,uj) sin out (7) 

The out-of phase response contains the information about ex¬ 
citation spectrum of the unperturbed system. In the next sec¬ 
tion, we will define the response function {k,iuj), where iuj 
is the Matsubara frequency, and relate this to the single parti¬ 
cle Green’s function of the bosons and hence to the excitation 
spectra. The imaginary part of IVa{k^ oj + then measures 
the amplitude of the out of phase modulation of the spin de¬ 
pendent momentum distribution due to the perturbation. The 
structure of the response is qualitatively different in the Mott 
and superfluid phase owing to the broken 17(1) symmetry in 
the superfluid phase and associated anomalous propagators. 
Hence we will treat the response in the Mott and superfluid 
phase separately. 


III. RESPONSE IN THE MOTT PHASE 


In the strongly interacting limit, the system is in an incom¬ 
pressible Mott insulating state with gapped excitation spec¬ 
trum. In the atomic limit, (7 = J = 0), the ground state has 
a fixed number of bosons of each spin, at all the sites. In 
general, is determined by the values of p and O. We will 
restrict our analysis to the case where n° = 1 and = 0, but 
the analysis can be easily extended to arbitrary integer values 
of n°. The single particle Green’s function in the Mott phase 
is a 2 X 2 matrix in the spin space, in terms of which the re¬ 
sponse function is given by 

IVa{k, iuJn) = ^ '^[G{k, iuJi)A{k)G{k, iwi + iuJn)]aa 

+ (^Ura —^n) ( 8 ) 


The boson Green’s function can be worked out in a strong 
coupling expansion around the localized atomic limiP^ElSIl 
In the Mott phase, it is given by 


G '^{k,ioJn) 


Fi{i 


^n) 

-l*k 


-Ik \ 
F2iiuJn) - ^k ) 


Fi[iu>n) — iuJn Eq — 2U -f 2C /^/(ilOn F Eq U), 
F 2 {iuj„) = ioJn - El (9) 


where Eq = fi + fl and Ei — ft + C,U — li¬ 
lt is evident that this Green’s function is not diagonal in the 
basis of non-interacting bands, as the atomic limit local prop¬ 
agator is different for the 2 spin species. This can be traced 
to the fact that in the atomic limit (J = 0, 7 = 0) the H term 
lifts the degeneracy between the spin states and as a result one 
obtains a polarized Mott state which is no = 1 for the f spins 
and no = 0 (vacuum) for the ), spins. The Green’s function 
can, however, be diagonalized to obtain 


GoikGuJn) 


0 



C+(k,iuj„) J 


( 10 ) 
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FIG. 5. Excitation spectrum in the superfluid phase. The parameters 
of the system are J = 0.03t/, 7 = O.OIG, p = 0.351/, H = 
0.011/ and = 0.4. In this case the BEC is formed at the zone 
center. The excitations are (a) two phase modes (b) one phase and 
one amplitude mode (c) one phase mode and one interband transition 
(d) two amplitude modes and (e) an amplitude mode and an interband 
transition. Note that the color scheme has same absolute value in all 
panels to show which excitations overlap with each other in energy. 


where the diagonal components are given by 

C,±{kGuJn) = E+{iLkln) - efc ± iEl{iu}n) + \lk\‘^y^‘^, (11) 


and F±{iijJn) = (l/2)[Fi(/u;„) ± F 2 {iujn)]- The basis trans¬ 
form, which diagonalizes the Green’s function is given by 


y-v {k,lU}„) U {k,lUJn) J 

|■u(A:,^w„)p = ^l/N{k,iu}n) 

v{k,iuJn) = 'yj:^u*{k,ioin)[Fi{iuJn) - - C+{k,iuJn)], 

N{k,iu}n) = iJ[Ei{iuin) - efc - C-i-(fc, *w„)]2 - 1 - 7 ^. (12) 


Note that -f |u(fc,la;„)p = 1. 

It is evident from the above expressions that even in the 
transformed basis, the diagonal Green’s function (Eq. [T0| ) 
will have a complicated frequency dependence. However, we 
are only interested in the out of phase response of the sys¬ 
tem to the external perturbation, which depends solely on the 
imaginary part of the Green’s function, analytically contin¬ 
ued to the real frequency domain (iw„ —w + 10+). It can 
be easily shown that the retarded Green’s function only has 
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simple poles, and so, for the purpose of calculating the out of 
phase response, the complicated expression of Eq. [T^can be 
replaced by the simpler form 


Goik^iun) 


iUJn—E~^ 


V 


0 


0 


z 




iUJn—E 


+ 


z 


+ 

hk 




(13) 


where uj = are the two zeroes of Q'^{k,uj) and corre¬ 

spond to the particle and hole excitations of the -f band, while 
oj — E~f. is the zero of (~{k,uj). The residues z are given 
by the oj derivative of the components of the Green’s func¬ 
tion, Eq. 10 evaluated at the corresponding poles. The three 
poles have a relatively simple explanation in terms of the par¬ 


ticle and hole excitations of a Mott insulator. In absence of 
the spin-orbit coupling (7 = 0 ), the f spins form a no = 1 
Mott insulator and has the corresponding particle and hole ex¬ 


citations. This corresponds to the limiting form of E 


p{h)k 


as 

7 —> 0 . The 4 , spins , however, form a no = 0 Mott insulator, 
and hence has only particle excitations. This is the limiting 
form of E~f, in the limit 7 —>^ 0. In presence of spin-orbit 
coupling, the spin states get mixed, but this three pole struc¬ 
ture of the Green’s function persists. We note that the Green’s 
function in Eq. 13 has the same imaginary part as the actual 
Green’s function (Eq. [T0| and hence there is no approxima¬ 
tion made in the above replacement, as far as computation of 
the out of phase response of the system is concerned. 

The response of the system can now be calculated in terms 
of the diagonal Green’s function as 


Ila{k,iU}n) 


A 

'p 


^M{k,wi)GD{k,i^i)M ^{k,ioJi)K.{k)M{k,ioJi + ioJn)GD{k,ioJi + ioJn)M ^(fc, iu;;-f 

ioJi 


GG 


EoJn —>• —UJn, 


(14) 


where A(k) is given by Eq. Note that the transformation 
matrices are themselves function of the Matsubara frequen¬ 
cies icji and iun- However, it can be easily seen that the trans¬ 
formation matrix elements, analytically continued to the real 
frequency domain {iuji —> w -f iO'*'), have no imaginary part. 


and hence the real frequency response will be dominated by 
the singularities of the diagonal Green’s function only. 

Analytically continuing to real frequencies and working out 
the Matsubara sums (and noting that matrix elements of M are 
real in this limit), we get 


,(2)t 


/ OO 7 / 

— apg{k,oj',uj + oj')G^jj {k,oj')G]j {k,uj + uj')[nB{uj') - UBioo' + w)] (15) 

-00 E 

where p, g = ±, G^ indicates the imaginary part of the Green’s function component, and the matrix element 


apqik,0Ji,0J2) = ^M^p{k,0Ji)Mp^{k,UJi)Amnik)Mnq{k,UJ2)MgJ{k,U}2) 


(16) 


We now restrict ourselves to the response at T = 0. Erom the simple pole structure of the Green’s function the response is then 
obtained as 


>i(7,k,u) = 27rA [a^+ik, E^^^, - E^^ + E^^^) + al_{k,E^^,E^^)zt^f^Zpf^5{oj - E^^ + E^,^) 

+Q^++(^) Ehk^ ~ ^hk + ^vk) + ^^hk’ ^nk)^hk^r>k^i^ ~ ^hk + ^vk) 


(17) 


It is clear from Eq.[^that the response shows up on two con¬ 
tours specified by w = E+^ - E+^ and 0 : = E-^ - £;+, 
the intra and interband particle-hole excitations. In the Mott 
phase, both these excitations are gapped with the lowest Mott 
gap corresponding to the intra-band excitations. As the fre¬ 
quency of modulation is swept, there is no response till the 
frequency crosses the Mott gap. Beyond this point, two dis¬ 
tinctive set of phenomena can be seen depending on the ratio, 
7 /J. Eor small 7 /J the minimum of the intraband excitations 
is at the zone center, [0, 0]. This is seen in Eig. [^a), where 


the intraband excitation in the Brillouin zone is plotted as a 
color plot for the following set of parameters: J = 0.02(7, 
7 = 0.01(7, n = 0.01(7, p = 0.2(7 and C = 0.4, where the 
minimum of the spectrum (the Mott gap) is 0.24(7 at the zone 
center. The dispersion increases as one moves away from the 
zone center and the highest excitation energy occur at [tt, tt 
with OJ = 0.57(7. The corresponding inter-band excitations 
are shown in Eig. [^b). They follow a similar pattern as the 
intra-band excitations with a minimum energy of 0.73(7 and a 
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maximum energy of 1.22U. 

The optical modulation response in this case is plotted for 
four different frequencies in this case in Fig. As frequency 
increases beyond the gap, 0.24[/, a contour of excitations is 
seen, which become larger and moves towards the edge of 
the Brillouin zone, before disappearing at w = 0.5717. This 
contour, corresponding to intra-band excitations, are shown in 
Fig|^a) and (b). As the frequency is further increased beyond 
u! = 0.73C/, the contours corresponding to the inter-band ex¬ 
citations appear around the zone center and move outwards, 
before disappearing at w = 1.2217, as shown in Figj^c) and 
(d). We would like to note here that the full dispersion of the 
excitations can be tracked from the optical lattice modulation 
spectroscopy. To show how this is done, in Fig. we plot the 
lattice modulation response in the Mott phase as a function of 
frequency along various cuts in the Brillouin zone, going radi¬ 
ally outwards at different angles 9 with the axis. Fig. ref- 
figS (a)-(c) correspond to the system parameters J = 0.0217, 
7 = 0.0117, = 0.217,n = 0.0117 and ( = 0.4 with the band 

minimum at the zone center, which is clearly seen from the 
plots. Further, the modulation spectroscopy also provide in¬ 
formation about the spectral weight of the various excitations. 
In this case, it is evident the interband excitations carry much 
less spectral weight than the intraband excitations. 

The phenomenology changes dramatically if the spin or¬ 
bit coupling 7 is larger than hopping J. Fig. [^c) and (d) 
shows a plot of the intra and interband particle hole excita¬ 
tions in the Mott phase for a system with parameter values 
J = O.Olf/, 7 = 0.0417, n = 0.0117, /i = 0.21/ and C = 0.4. 
In this case the spin orbit coupling is larger than the hopping 
and consequently the intraband excitations have a minimum 
at [±fco)ifco] along the zone diagonal. The location of the 
minima corresponding [Aq, ko] is shown in Fig. [^c); the cor¬ 
responding figures for other minima can be obtained by a re¬ 
flection of the contours in the first quadrant (0 < kx,ky < tt) 
about the k^ and ky axes and the origin. The interband exci¬ 
tations however have a minimum at the zone center, as seen 
in Fig. [TJd). Once again, the system will not show any re¬ 
sponse as long as the modulation frequency is below the Mott 
gap of 0.281/. Once the Mott gap is crossed, excitation con¬ 
tours around appear in the response, as seen (for 

[ko,ko]) in Fig|^a). These contours spread out and disap¬ 
pear at ct) = 0.451/, as seen in Fig [^b). Another contour, 
corresponding to interband excitations and centered around 
k = 0 appears as the modulation frequency is swept beyond 
0.871/. This spreads out and finally disappears atuj = 1.11/, 
the upper limit of the interband excitation energy, as seen in 
Fig. I^c) and (d). To get a better idea of the dispersion, the 
lattice modulation response for a system in Mott phase with a 
band minimum shifted to finite wavevectors is plotted in Fig. 
El(d) - (f), where the parameters of the system are J = 0.011/, 
7 = 0.041/, /i = 0.21/,n = 0.011/ and C = 0.4. In this case, 
the band minimum is at {kQ,ko) with fcg = 1.131a“^. It is 
clearly seen that the intraband excitations have a shifted band 
minimum, while the interband excitations continue to have a 
minima at the zone center. In this case, it is clear from compar¬ 
ison with corresponding response in Fig. that although the 
intraband excitations have larger spectral weight than the in¬ 


terband excitations, the interband excitations carry relatively 
larger spectral weight than the case where the minimum was 
at the zone center. 

The lattice modulation spectroscopy can identify the Mott 
phase from the existence of a gap in the spectrum. How¬ 
ever the momentum resolved nature of the spectroscopy pro¬ 
vides much more detailed information about the single parti¬ 
cle spectral function of the system. It provides information 
about the spectrum and one should be clearly able to see the 
minimum of the excitation spectrum shift from the zone cen¬ 
ter as the parameters are varied. The modulation spectroscopy 
also provides information about the spectral weight of the ex¬ 
citations of the system, which is of immediate relevance for 
figuring out both near equilibrium and far from equilibrium 
response of the system to different stimuli. 


IV. RESPONSE IN THE SUPERELUID PHASE 

The spin-orbit coupled bosons undergo a Mott insulator- 
superfluid quantum phase transition as either the hopping or 
the spin-orbit coupling is increased as both terms help to de¬ 
localize the bosons. If the transition occurs at a large value 
of J/ 7 , the superfluid phase has a BEC at the zone center 
k = [0,0] with associated anomalous propagators and Gold- 
stone modes. On the other hand, if the transition takes place at 
a large value of 7 /J, the system exhibits the twisted superfluid 
phase, with a BEC at a finite momentum. In general there can 
be four such momentum values given by [±7o, ±^o] which are 
degenerate minima of the spectrum. This in principle allows 
the possibility of formation of a square superlattice (i.e. a su¬ 
persolid) with the incommensurate momentum ko providing 
the inverse lattice constant. However, in a cold atom system 
with the presence of a trapping potential which breaks trans¬ 
lational symmetry, the more likely effect is the formation of 
domain^, each of which corresponds to choosing one of the 
four values of allowed ko. Assuming domains of size much 
larger than k^^, the momentum distribution signal will be an 
incoherent weighted sum of the signal from a single domain 
with a fixed condensation wavevector. In this paper, we shall 
assume that condensation occurs at only one possible momen¬ 
tum which we choose to be {ko, ko) and work out the signal 
from a single domain. The important features, like the low 
energy features around (7o, ^o) will occur in different parts of 
the Brillouin zone for different domains and thus will not be 
washed out by formation of domains. Although the detailed 
signal requires knowledge of domain distribution, a first ap¬ 
proximation can be obtained from our calculation by symme¬ 
try, together with an idea of the density of each type of do¬ 
mains. The lattice modulation spectroscopy is uniquely suited 
to distinguish between gapless superfluid phase from the Mott 
phase, and owing to its momentum resolved nature, it can also 
distinguish between the normal and twisted superfluid phases. 
Since this method gives detailed information about the spec¬ 
tral function, useful quantities like the speed of sound can be 
calculated from the spectrum, while the relative intensity of 
the contours will provide information about transfer of spec¬ 
tral weight among the different kinds of excitations as the en- 
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FIG. 6 . Lattice modulation response in the superfluid phase. The parameters of the system are J = 0.0317, 7 = O.Olf/, /r = 0.35C/,G = 
O.Olf/ and ( = 0.4. In this case the EEC is formed at the zone center, (a) A single contour corresponding to the two phase modes (b) 2 
contours (outward from zone center): a phase and an amplitude mode, two phase modes (c) 3 contours (outward from zone center): a phase 
mode and an interband transition, a phase and an amplitude mode, two phase modes (d) 2 contours (outward from zone center): a phase mode 
and an interband transition, a phase and an amplitude mode (e) 3 contours (outward from zone center): two amplitude modes, a phase mode 
and an interband transition, a phase and an amplitude mode (f) 2 contours (outward from zone center): two amplitude modes, a phase mode 
and an interband transition (g) 3 contours (outward from zone center): an amplitude mode and an interband transition, two amplitude modes, 
a phase and an interband transition (h) 2 contours (outward from zone center): an amplitude mode and an interband transition, two amplitude 
modes (i) a single contour: an amplitude mode and an interband transition. Note that the spectral weight of the interband transitions are much 
smaller than the weights of other transitions. 


ergy is varied. 

In the superfluid phase the bosons occupying the lower en¬ 
ergy band form a EEC at the appropriate momentum. In this 
case the Green’s functions expand to a 4 x 4 matrix to ac¬ 
commodate the anomalous propagators. It is easier to work 


in the ± basis, earlier described in the Mott phase, since the 
structure of the Green’s functions are simplest in this basis. 
Working with a 4 component vector [(j>~ {k, iui), cj)~* {2kQ — 
k, —iuj), (l)~^{k, iuj) .(l)^*(2 kci — k, —iuj)], the inverse Green’s 
function is given by^oEIIIl 




C-{k,iUJn) 

0 

0 

0 


0 0 

(2fco — k, —iuJn) — r 0 

0 C+(fc,zw„)-r 

0 -r C-r(2^o 


0 

0 

—r 

k, -iuJn) 



( 18 ) 


where (A:, iuj) is defined by Eq. 11 and r = C+(fcoj ko,uj = 
0 ) incorporates the effects of the presence of the condensate. 
The lower 2x2 block simply looks like the Green’s func¬ 


tion of a Bose gas in the Bogoliubov approximation, with the 
complicated function C-i-(^) incorporating correlations due 
to proximity to a Mott insulator, replacing a simple free par- 
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tide piece. In the diagonal upper block, the presence of the 
condensate leads to a Hartree correction due to inter-band in¬ 
teractions. 


FIG. 7. Excitation spectrum in the superfluid phase. The parameters 
of the system are J = O.Olf/, 7 = 0.06(7, p = 0.2(7,17 = 0.01(7 
and ^ = 0.4. In this case the EEC is formed at a finite momentum 
[±A:o,±fco], with fco = 1.34a“^, away from the zone center. The 
excitations are (a) two phase modes (h) one phase and one amplitude 
mode (c) one phase mode and one interband transition (d) two ampli¬ 
tude modes and (e) an amplitude mode and an interband transition. 
Note that the color scheme has same absolute value in all panels to 
show which excitations overlap with each other in energy. 


The unitary transform which converts the original “spin” 
basis, i.e. ioj), — k, —iuj), </>^(fc, ioj), ^f>|(2fco — 

k, —iuj)], to the ± basis is given by 


M{k,iujn) 


/ Mil(kjiuj) 0 M2i{k,iuj) 0 

0 Mii{2ko — k,—iw) 0 M2i{2ko — k,—iuj) 
Mi2{k,iuj) 0 M22{k,iuj) 0 

\ 0 Mi2{2kQ — k,—ioj) 0 M22{2ko — k,—iuj) 


where the matrix elements Mij are given by Eq.[T^ The final Green’s function has the form 
/ ( 

G{k,iuJn) = 


G-{k,iojn) 0 0 0 

0 G~ {2kQ — k,—iujn) 0 0 

0 0 G+{k,ii^n]ko) F{k,iuJn]ko) 

0 0 F(2fco — fc, —fco) G'^{2kQ — k,—iuJn',kQ) 


(19) 


( 20 ) 


The full expression for the Green’s functions are complicated, of the Green’s functions are simple poles. So, as far as the 

but, as in the case of the Mott phase, the only singularities imaginary part of the Green’s function is concerned, they are 

faithfully reproduced by 
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kx 


FIG. 8 . Lattice modulation response in the superfluid phase. The parameters of the system are J = O.Olf/, 7 = 0.0617, /r = 0.2I7,fI = 0.0117 
and ( = 0.4. In this case the BEC is formed at a finite momentum [fco, ^o] (fco = 1.34a“^) away from the zone center. A single contour 
corresponding to the two phase modes (b) 2 contours (outward from zone center): a phase and an amplitude mode, two phase modes (c) 4 
contours : two open contours due to two phase modes, and two closed contours (outward from the zone center) due to two amplitude modes, a 
phase and an amplitude mode (d) closed contour due to a phase and an amplitude mode and open contours due to two amplitude modes, (e) a 
single contour: two amplitude modes (f) a single contour: a phase mode and an interband transition at finite momentum (g) 2 contours (outward 
from zone center): an amplitude mode and an interband transition, a phase and an interband transition (h) a single contour: an amplitude mode 
and an interband transition. Note that the spectral weight of the interband transitions are much smaller than the weights of other transitions. 


G {k,ioj) 


F{k,iu}-,ko) 


9pk 


- Epk ’ 

fik 

iu - 




G+{k,tu;;ko)= . + 


Jlk 

iu) + 


+ 


J2k 

iu: - E+^ 


atk 


lu: - E+^ 


lUJ - 


Hk 


■Kk 


+ 




iu: - E+^ 


iu: + E+^^ 


atk 


iu + E+^ 


( 21 ) 

( 22 ) 


Here the particle excitation in the — branch, 77^^ is the solu¬ 
tion of C-(fc) Epk) — '^he four excitation poles ±i7^2)fc 
are obtained from the solutions of 

[C+(/c,a;) - 2r][C+(2fco - fc,-w) - 2r] - = 0. (23) 

The quasiparticle residues are obtained from the frequency 
derivatives of the Green’s functions at the corresponding 
poles. Note that in the + band, the particle and hole excita¬ 
tions are now mixed due to formation of a condensate and that 


the energies T^^ 2 )fc’ 5 i( 2 )fc^’ ft( 2 )k depends on /cq through 
Eq.|^ this is in contrast to the corresponding quantities in the 
— band which has no particle-hole mixing. This necessitates 
the use of additional fco argument in the dehnition of G~^ and 
E in Eq.|^ however, we refrain from putting such additional 
label of fco in £^^ 2 )fc’ 5142 )^^’ -^ 1 ( 2 )^ notational brevity. 

The response function for the lattice modulation spec¬ 
troscopy is then given by 


11(7 (fc; i^n) 


A 


|^A^(fc, iu:i)G{k, iu:i)M ^(fc, iu:i)A {k)M{k, iuii -\- iu:n)G{k, iui + iu:n)M ^(fc, iui + iun) 

ibJi 


2a—1,2a—1 


EUn —Un 


(24) 
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(a) 0=71/2 (b) 0=71/4 



(c) 0=7r/6 



FIG. 9. Lattice modulation response in the superfluid phase with the variation in frequency for various cuts in the Brillouin zone going radially 
outwards at an angle 6 with the k^. axis. In (a) - (c) the parameters of the system are J = 0.03?7, 7 = 0.0117, /i = 0.3517,17 = 0.0117 and 
= 0.4 with the condensate at the zone center, (a) 9 = 7 r/ 2 , (b) 9 = 7 r /4 and (c) 9 = -k!%. All the five branches including the gapless linear 
Goldstone mode near fc = 0 is clearly seen. In (d) - (f) the parameters of the system are J = 0.0117, 7 = 0.0617, = 0.217,17 = 0.0117 and 

= 0.4. In this case the condensate is located at (fco, Ico) with fco = 1.131a“^. (d) 9 — 7 r/ 2 , (e) 9 — 7 r /4 and (f) 0 = 7 r/ 6 . Only the cut 
along [1,1] direction passes through the condensate location and shows the gapless Goldstone mode, while the spectrum is gapped along the 
other cuts. 


where a = 1 for cr =t and a = 2 for a = 4 ,, and the perturbation matrix A' is given by 


A'(fc) 



/ efc 

0 

7fc 

0 

1 

0 

efc 

0 

* 

7-fc 

2 

7fc 

0 

efc 

0 


V 0 

7-fc 

0 

efc 


(25) 


Working out the Matsubara sums, and analytically continuing to real frequencies, 


fc,a;) = 2A 


/ oo 7 / 

—w + {k,uj'; ko)G‘^ {k,uj + uj'; ko)[nB{uj') - nsioj'+ uj)] 

-00 TT 


+P'^ {k,uj',uj + u!')F {k,u;';ko)F (fc, w — w'; fco)[ti-B(w') — — w') 


(26) 


where p,q = ±, and F indicate the imaginary part of the Green’s function component, and it is understood that 
G~ {k, oj; fco) = G~ (fc, uj). Here the matrix elements 

apg{k,uji,uj2) = '^M2a-ipik,uJi)M~^{k,uJi)A,^„{k)Mnq{k,uj2)M~2a-i{k,uj2) 

mn 

l3‘^{k,UJi,UJ2) = ^M2<,-l,3(fc,U;i)M2;/_4(fc,Wi)A2„_2m((':)('/^2n.4(fc,W2)M^2'cr-1.3(*>‘^2) (27) 

mn 

We now restrict ourselves to the response at T = 0 for a; > 0. The w < 0 response can then be obtained from the fact that 
the imaginary part of the response function is an odd function of uj. From the simple pole structure of the Green’s function the 
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response is then obtained as 


,( 2 ) 


{a,k,uj) = 2 tt\ \^p‘l{k)6{uj - - E+^) + p^{k)S{uj - - E+^) 

pl{k)5{uj - 2 ^;+) + p'l(k)5{u: - 2E+^) + p'^{k)S{uj - £+ - E+^)] 


( 28 ) 


where the weight of the different contours are given by 


PlW — ‘^9pk9lk^l2i^T^pk^ ^Ik) 

pm = 2g+,^gt,mmm-Etk) 
pm = 2g+^g+^aUk,E+^-E+k) 
pm = 2gt^gtk^aUk,Et„-E+k) 


P2i^) — ‘29pk92k^l2i^yEpi^, E^i.) 

(fm[l3^ik,Et„Et,)+l3^ik,mk^-Etk)] (29) 

(/+)2r (A, E+„E+,) + ^%k,-E^, -E+k)] 

‘^9ti9tka^22{K - 2/i+,/2+,r (fc, E^,,E^,) + ^^k, -E+k^ -^tk)] 


It is clearly seen that the system would show response on the 
contours which correspond to energies 2E^f,, 2i?^,, + 

E^k, E~f. + E^j^ and E~f, + E^i,. Here E^^, is the Goldstone 
phase mode which goes gapless, is the gapped amplitude 
or Higgs mode and E~f, is the particle excitation to the up¬ 
per band, which also has an excitation gap. The location of 
these contours in the Brillouin zone are dramatically different 
depending on whether the BEC is formed at the zone center 
(large J/7) or at a finite momentum (small J/7). 

The dispersions of 2E;+, 2E^f,, + £;+, E~^. + and 

E~k + E^i. are plotted in Fig. ^for a system with following 
set of parameters: J — 0.03(7, 7 = 0.0H7, Q — 0.01(7, p — 
0.35(7 and = 0.4. In this case, the system is in a superfluid 
with a BEC at the zone center. Fig. |^(a) shows 2i?^, which 
ranges from 0 at the zone center (the gapless point) to 0.58(7. 
All the other plots in Fig. show similar features except the 
fact that these excitations are gapped. In Fig. §(b), E++E+ 
ranges from 0.37(7 to 0.83(7. Similarly E . + E^^ [Fig.l^ 
(c)] ranges from 0.51(7 to 1.2(7, 27(7 [Fig. [^(d)] ranges from 
0.74(7 to 1.1(7, and E~^-\-E^f. [Fig.^(e)] ranges from 0.88(7 
to 1.5(7. ^ 

The optical modulation response for the above system is 
shown in Fig|^ Fig. |^a) shows the response at w = 0.2(7, 
where a single contour of excitations correspondi^ to excit¬ 
ing two phase modes, uj — 27(^ is seen. Fig. ^b) shows 
the response at a; = 0.4(7, where two contours of^xcitations 
are seen, the outer one corresponding to exciting two phase 
modes, w = 2 £'^, and the inner one corresponding to excit¬ 
ing a phase and an amplitude mode uj = E^f. + E^f^. Fig.^c) 
shows the response at w = 0.55(7, where three contours ofex- 
citations are seen, the outermost one corresponding to exciting 
two phase modes, oj = 2 £(^, the middle one corresponding 
to exciting a phase and an amplitude mode w = 7(^ + E^f., 
and the innermost contour corresponding to exciting a phase 
mode and a band transition to the — band, to = E^kEEpf.. As 
the frequency is increased to w = 0.6(7, the two phase modes 
disappear and two contours corresponding to w = E^k + E^f, 
and OJ = E^f. -f E~f, remain, as seen in Fig. |^d). As the fre¬ 
quency is further increased to w = 0.8(7, an additional con¬ 
tour (the innermost) corresponding to w = 27([^ appears in 
Fig.|^e) along with the other excitations seen in Fig. |^d). At 
UJ = 0.85(7, the contours in Fig|^f) corresponds to w = 2E^j^ 


and OJ = Etk + Epk, while at w = (7, Fig|^g) has an ad¬ 
ditional innermost contour of w = + 7(“^. Finally, in 

Fig|^h), there are two contours corresponding to w = 2E^j, 
and OJ = E^i, -f E~f^ at oj = 1.15(7, while at w = 1.4(7 
(Fig.^i) ), only the contour corresponding to w = E^^ + E~f^ 
remains.The presence the gapless phase mode as well as four 
other gapped modes is clearly seen in Fig. (a)-(c), where 
the lattice modulation response is plotted as a function of fre¬ 
quency along different radial cuts in the Brillouin zone mak¬ 
ing angles 6 = 7r/2 , 0 = 7r/4 and 6 = tt/6 with the kx 
axis respectively. The speed of sound in the system can be 
obtained from the slope of linearly dispersing phase modes 
seen in the figure. The gapless mode and the unique pattern 
of obtaining up to three contours at certain frequencies distin¬ 
guishes the superfluid phase from the Mott phase and can be 
used to detect the superfluid-insulator quantum phase transi¬ 
tion in this system. One can also obtain detailed information 
about the spectrum and relative spectral weights of the various 
modes from the lattice modulation response. 

We now consider the changes that appear in the opti¬ 
cal modulation spectroscopy when the spin-orbit coupling is 
greater than the hopping amplitude and a BEC is formed at a 
finite momentum. To see this, in Fig. |7] we plot the disper¬ 
sions of 2E^f., 2i?^, E^i^ + E^f,, Epi^ + E^f, and E^j^ + E^j^ 
for a system with following set of parameters: J = 0.01(7, 
7 = 0.06(7, H = 0.01(7, p = 0.2(7 and C = 0.4. In this 
case, the system is in a superfluid with a BEC at the momen¬ 
tum [1.34a“^, 1.34a“^] where a is the lattice constant. Fig.[^ 
(a) shows 27^^, which ranges from 0 at BEC wavevector (the 
gapless point) to 0.4(7. All the other plots in Fig.|7]show sim¬ 
ilar features except the fact that these excitations are gapped. 
In Fig.M(b), 7 ;+ -f 7 ;+ ranges from 0.1(7 to 0.42(7. 27;+ 
ranges from 0.21(7 to 0.51(7 as shown in Fig. I^(c). The ex¬ 
citations involving band transitions, E~f^ -\- E^ ranges from 
0.86(7 to 1.06(7 [Fig.|^(d)], and 7;"^ +7;+ [Fig.|^(e)] ranges 
from 0.9(7 to 1.13(7. 

The optical modulation response for the above system is 
shown in Fig[^ Fig.j^a) shows the response at w = 0.02(7, 
where a single contour of excitations corresponding to excit¬ 
ing two phase modes, oj = 27;^ is seen around the BEC 
wavevector [ko,ko]. Fig. [^b) shows the response at w = 
0.15(7, where two contours of excitations are seen, the outer 
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one corresponding to exciting two phase modes, w = 
and the inner one corresponding to excit ing a phase and an 
amplitude mode w = Fig. ^c) shows the re¬ 

sponse at w = 0.2517, where three contours of excitations are 
seen, the outermost one corresponding to exciting two phase 
modes, w = 271^, the middle one corresponding to excit¬ 
ing a phase and an amplitude mode w = E^/. + and 
the innermost contour corresponding to exciting two ampli¬ 
tude modes, ui — 271^. As the frequency is increased to 
uj = 0.4217, the two phase modes disappear and two contours 
corresponding to a; = -I- £1^ and oj = 2i?^ remain, as 

seen in Fig. [^d). As the frequency is further increased to 
w = 0.517, only the contour corresponding to w = 2iil^ ap¬ 
pears in Fig. [^e). As modulation frequency is increased fur¬ 
ther, the system shows no response, till the frequency reaches 
w = 0.8617. Beyond this point, a contour of excitations corre¬ 
sponding to a; = E^f^+E~f^ appears around the BEC wavevec- 
tor, as seen in Fig.|^f) for to = 0.8717. Beyond to = 0.917, a 
second contour, w = E^^ + E~f. appears in Fig.j^g). Finally 
Fig.j^h) shows the response at w = 1.117, where a single con¬ 
tour corresponding to w = E^^^ + £1“^ is present. The distinct 
pattern of the spectrum is better visualized in Fig. (d)-(f) 
where the lattice modulation response is plotted as a function 
of frequency along different radial cuts in the Brillouin zone 
making angles 0 = 7r/2,0 = 7r/4 and 0 = tt/G with the 
axis respectively. For the 9 = tt/2 and 9 = nj6 cut, which 
does not pass through the condensate wavevector, all the spec¬ 
tra are gapped, while the 0 = 7r/4 cut, which passes through 
the condensate location shows the gapless Goldstone mode. 
This clear qualitative distinction along various cuts can lead 
to a precise location of the condensate wavevector. In addi¬ 
tion information about spectral weights can also be obtained 
from the lattice modulation signal. 


V. DISCUSSION 

In this work, we have studied a system of two species of 
interacting bosons in an optical lattice to modulation of the 
lattice potential. The spin states of the bosons are coupled 
through a spin-orbit coupling, which is implemented either 
by Raman dressing or by time-dependent magnetic field gra¬ 
dients. This system shows a superfluid-Mott insulator quan¬ 
tum phase transition as a function of increasing interaction 
strength. In addition, as a function of the relative strength of 
the spin-orbit coupling to the hopping amplitude, the system 
in the superfluid phase shows a transition from an ordinary 
superfluid phase with a BEC at the zone center to a twisted 
superfluid phase with a BEC at a finite momentum. We have 
provided a technique for differentiating between the different 
phases of such a system based on the response of such a sys¬ 
tem on to a modulating optical lattice. 

In addition to finding the location of the precursor peaks in 
the MI phase and condensate position in the SE phase which 
can also be obtained by other standard experimental tech¬ 
niques such as time-of-flight measurements. However, lat¬ 
tice modulation spectroscopy provides several additional in¬ 
formation. Eirst, one can use this technique to map out the 


single-particle excitation spectrum of the bosons. In the Mott 
phase, it provides the effective mass of the dispersion around 
the minimum gap. In the superfluid phase, this technique 
not only shows the presence of gapless excitations, a quan¬ 
titative estimate of the speed of sound Vs (the slope of the 
gapless mode) and the mass of the Higgs (gapped) mode can 
be extracted from the experimental data. Eurther, since the 
matrix elements of the lattice modulation operator can be ex¬ 
plicitly calculated within our technique, the spectral weight of 
the different modes can also be extracted from the modula¬ 
tion spectroscopy response. In this context, we would like to 
note that the expression of the response in terms of the spec¬ 
tral function is more general than the particular approximation 
used to calculate it in this paper. Eor example, in general the 
Higgs mode will develop a width due to decay to two phase 
modes. This can also be computed from the modulation spec¬ 
troscopy response. Einally, we would like to point out that 
there is a technical advantage of our calculation over its coun¬ 
terpart for Bragg spectroscopy. This advantage stems from 
the fact that the optical modulation response, being a zero mo¬ 
mentum transfer process, does not receive large contribution 
from the vertex correction terms. Thus, even if the poles in 
the single particle Greens function are broadened due to self 
energy corrections, the optical modulation response function 
will pick out the frequency convolution of the one particle 
spectral function at the same momentum. Thus this specific 
spectroscopic method provides direct access to single particle 
excitations in the system. 

The experimental verification of our theory invo lves use of 
standard spectroscopy experiment technique^^HH] 'pjjg ^pg. 
ciflc experiment that we propose involves modulating J and 
7 by a laser creating an additional optical lattice with mod¬ 
ulation frequency. After this modulation, one turns off the 
trap and the lattice and measures the position distribution of 
the outgoing bosons as done in any standard time-of-flight 
measurement. The position distribution of these bosons un¬ 
der standard experimental condition^ reproduces their mo¬ 
mentum distribution ni„od(k, w) inside the trap (in the pres¬ 
ence of the modulating lattice). We suggest a comparison of 
'^mod (k , w) to the momentum distribution n(k) of the bosons 
without the modulating lattice potential to obtain 5 n(k, w) = 
?T-mod(k, w) — n(k). We expect ( 5 n(k, w) to provide the neces¬ 
sary information about the boson spectral function and carry 
the signature of the Mott and the superfluid states as described 
in Secs. m and We note that realization of these ex¬ 
periments requires that the thermal smearing of the contours 
would be small enough to distinguish between the different 
phases. In the Mott phases this requires ksT 0.2(7, where 
ks is the Boltzman constant and ksT* = 0.2(7 is the melt¬ 
ing temperature of the MI phasd^. This is readily achieved in 
standard experiments where (7 ~ 2 —5 KHz ~ 200—500 nK. 
In the SE phase, the precise nature of the low-energy con¬ 
tours for the Goldstone modes may be difficult to discern since 
it would require a small T. However, the amplitude modes 
which occurs at finite energy scale would still be observable 
in a straightforward manner within current experimental reso¬ 
lution. An estimate of the maximal allowed thermal smearing 
for the SE phase comes from the criteria T <C where 
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is the critical temperature of the SF phase which can be esti¬ 
mated to be ksTc — 2 :oMax[J, 7 ] ~ O.IU ~ 20 — 50 nK. We 
note that T ^ 1 nK have been achieved in standard ultracold 
atom experiment^. 

In conclusion, we have shown, via explicit computation of 
(5n(k, uj) in the strong coupling regime, that the response of a 
spin-orbit coupled Bose system to a modulated optical lattice 
can differentiate between (a) Mott and superfluid phases and 
(b) systems with finite momentum BEC and zero momentum 
BEC. The momentum resolved nature of the optical modula¬ 
tion spectroscopy, which provides information about the one 


particle spectral function of the system, resolves the super¬ 
fluid phase from the insulator phase by presence/absence of 
gapless Goldstone modes in the two phases. Eurther, the pat¬ 
tern of excitation contours appearing in the Brillouin zone as 
the frequency is tuned is distinct in the two phases and further 
helps to distinguish the phases. The momentum resolution can 
also resolve states with condensates at finite momentum from 
states with condensate at the zone center by looking at the 
location of the low energy excitations, which are always cen¬ 
tered around the momentum where the BEC forms. We have 
suggested concrete experiments which can verify our theory. 
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